Thermoelectric properties of the interacting two dimensional electron gas in the 

diffusion regime. 
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We demonstrate that kinetic coefficients related to thermoelectric properties of the two dimen- 
sional electron gas in the diffusive regime are strongly influenced by electron-electron interaction. 
As an example we consider the thermoelectric coefficients of the diluted two-dimensional electron 
gas in Si(100) MOSFET's in the presence of charged-impurity scattering. We find that the screening 
anomaly at q = 2kp, also responsible for Friedel oscillations, leads at low electron densities to a 
large change in the thermoelectric coefficient for the thermopower. 
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Since a long time thermoelectric effects in metals are 
considered in textbooks, see for instance the book by 
Ashcroft and Merminji Nevertheless, all known results 
of explicit calculations for two dimensional metals are 
made within the random-phase approximation and are, 
therefore, only valid for weak electron-electron interac- 
tion, which means for very high electron densities. 2 In 
a number of publications^ — the attempt was made to 
describe real experiments and these results were extrap- 
olated down to low electron density where the Coulomb 
interaction is strong. We argue in this paper that an 
important term for the calculation of thermoelectric co- 
efficients was neglected in this extrapolation. 

The present calculation demonstrates that the interac- 
tion leads to a strong change in the kinetic characteristics 
of the two dimensional electron gas, such that all text- 
book predictions for the Seebeck and the Peltier effect be- 
came invalid in the limit of low electron densities. From 
experiment^— and the theoryii - — it is well known that 
electron-electron interaction is the origin of a significant 
change in the electron dispersion law. The influence of 
interaction on kinetic coefficients is an additional effect 
which is not related to a change of the dispersion law of 
electrons. 

The presence of temperature T and electric potential 
gradients (electric field E) in the metallic system is ac- 
companied by the flow of an electric current I and a heat 
current J: 



I = L E 



3 = L 2 



L 12 (-VT), 
E + L 22 (-VT), 
L 21 = TL 12 . 



(1) 



In the heat transfer phonons prevail at high tem- 
peratures and thermoelectric effects are caused by the 
electron-phonon interaction. In the low temperature 
limit electrons give the main contribution to the heat 
flow and the thermoelectric effects at low temperatures 
are caused by the diffusion of electrons. 



Below we discuss thermoelectric effects due to diffu- 
sion in the two dimensional electron gas in the presence 
of electron-electron interaction. At low electron densi- 
ties n the interaction effects become very large. In the 
best two dimensional electron systems^ the value of di- 
mensionless Wigner-Seitz parameter r s = l/(7rna 2 3 ) 1 / 2 , 
describing the importance of the interaction effects, is 
about 10. Here as is the effective Bohr radius. 

Parameter r s — 10 means that the potential en- 
ergy exceeds the Fermi energy (kinetic energy) by 
about one order of magnitude. Recently very clear ex- 
pressed evidences of the strong electron-electron interac- 
tion were experimentally observed in Si-MOSFETs (see 
e.g. Refs^ - — ). Having this fact in mind we will apply 
our consideration to (100) Si-MOSFET with high mo- 
bility. Of course our results will be valid for other two 
dimensional electron and hole system. 

To the best of our knowledge thermopower calculations 
in the low density range with strong electron-electron in- 
teraction were never realized, despite the fact that the ex- 
perimental activity in two-dimensional electron systems 
is directed in lowering the electron density. We restrict 
our consideration to the region of low electron densities, 
where the system is still chacterized by a metallic con- 
ductivity a >> e 2 /h. We consider a simple model of a 
two-dimensional electron system in the metallic regime 
and we use the relaxation-time approximation. A gas 
of Landau quasiparticles has, renormalized by interac- 
tion effects, a mass m and a thermodynamic density of 
states In our calculation we neglect with heat trans- 
fer by low energy collective modest which has minor 
influence on the thermoelectric coefficients. Having in 
mind a low electron density system we suppose charged- 
impurity scattering by impurities randomly distributed 
at the Si/SiC>2 interface with an impurity concentration 
m. The extension effects of the electron wave function 
into the silicon are described by the formfactor Fi(q) 
for the electron-impurity interaction and by F(q) for the 
electron-electron interaction. 

The appearance of thermoelectric effects is caused by 
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the fact that the diffusion coefficient near the Fermi level 
is a function of the energy. A temperature gradient pro- 
duces two diffusion streams of electrons: "cold" electrons 
moves along the temperature gradient and "hot" ones in 
the opposite direction. If the diffusion coefficient is inde- 
pendent of energy, charge transfer is absent. In the case 
of an energy independent transport scattering time the 
diffusion coefficient increases with energy, which leads to 
a charge flow opposite to the temperature gradient, i.e., 
to the appearance of a charge flow typical for ideal elec- 
tron system. The contribution from an energy dependent 
transport scattering time should be added to this term. 

For a strong electron-electron interaction the screening 
function and the transport scattering time show a pecu- 
liar behavior at the Fermi energy. Scattering processes 
for electrons below and above the Fermi level are differ- 
ent due to the difference in the screening propertie o 17 i 18 
which can be also interpreted as additional scattering due 
to Friedel oscillations for high energy electrons^. On the 
first glance this peculiar behavior of the screening func- 
tion gives only small contributions to the transport scat- 
tering time and can be omitted for zero temperaturej 17 i 18 
However, the contribution to the energy dependence of 
the diffusion coefficient has the same order of magnitude 
as other terms and is therefore very important. This ef- 
fect was never considered in the past and will be shown 
in the following. 

Below we show that the " ideal electron system" contri- 
bution to the energy dependence of the diffusion coeffi- 
cient for a strong electron-electron interaction is compen- 
sated in part by the energy dependence of the transport 
scattering time. As a result the Seebeck effect and Peltier 
effect at low electron densities are expected to be signifi- 
cantly smaller than those for a non- interacting electron 
gas and even might change the sign^S 

Thermoelectric coefficients are sensitive to the energy 
dependence of the transport scattering time and are equal 
to 



Lf = ^r= 9 ^l I v fc v fe r (e) 
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2tt 



{e-e F ),dS u a/ , 



T 



\v k \ oe 
(2) 

where Vk is the electron velocity, dS — \k\d<fi, f° is the 
Fermi-Dirac distribution function, and g v is the valley 
degeneracy. 

The transport scattering time To(e) for a particle of 
energy e in lowest order of the random potential < 
\U(q)\ 2 > is given byi 7 -^ 
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To(e) 2nhe 



2k 



dq 



< \U(q)\ 2 > 



(4k 2 -q 2 )V 2 e(q) 2 



(3) 



where < \U{q)\ 2 >= mV (q) 2 F^q) 2 with V(q) = 
2ire 2 /eLq, and cl = \{^sc + tins) is the dielectric con- 
stant of the background material^ 

For a non-interacting electron gas with short-range 
random potential as disorder the time tq is independent 
of energy and only the " ideal electron system" contribu- 



tion to the thermoelectric coefficients exists. In this case 
the thermopower is given by the well known expression 



' do 
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L\2 
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T:g v k 2 B mT 
Sen 



k B T<^e F . (4) 



Here m is electron effective mass. 

With the aim to take into account the screening of 
impurities and especially the anomalous screening for q > 
2k F we write 



e(q) = e 1 (q)[l 



V(q)(l - G(q))p F F(q) 
61 (9) 

(i_^)V»e(«»-4^)], 



(5) 



where e(q) is the screeing function of the two dimensional 
electron system, G(q) describes the local-field correction 
for exchange and correlations effects not taken into ac- 
count in the random-phase approximation. 22 Q(x) is the 
step function with Q(x) = 1 for x > and Q(x) = for 
x < 0. ei(q) is defined by 



e 1 (q) = l + V(q)F(q)(l-G(q))p F 



(0) 



We use for the explicit calculation the local-field correc- 
tions in the Hubbard (H) approximation G(q) — Gn(q) 
where an analytical expression is available as 



G H (q) 



1 



2g v (q 2 + k 2 F y/ 2 - 



(7) 



For Si(100) MOSFET structures the valley degeneracy 
factor is g v = 2. 

With Eqs. (|5I6[) we can write inverse transport scatter- 
ing time as 
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Ms) T (e F ) ri(e)' 



(8) 



For low temperatures only terms proportional to (e — 
£f)/£f need to be taken into account in the energy de- 
velopment of the transport scattering time in Eq. |8] Only 
these terms will be important for the calculation of Lij 
for low temperatures. 

We evaluate the energy independent part as 



4m < \U(2k F )\ 2 > 



x dx 



r (£ F ) " nh 3 [ei(2M] 2 Jo [1 ~ x 2 ]V 2 

m < \U(2k F )\ 2 > 
~ ¥ [ei(2M] 2 ' 

The contributions to — consist of two terms: 
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m d < \U(2k)\ 2 > 
{£ ~ £F 'h 3 de [ [ £l (2fc)] 2 J 



6(e - e F )V(2k F )F(2k F )(l - G{2k F ))p F 



(e - e F ) 2m 

+ e F h 3 

< \U(2k F )\ 2 > 



el(2k F ) 



(10) 
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As we will see, the first is the regular term of the expan- 
sion and gives only a small contribution to thermoelectric 
coefficients at low electron densities. More important is 
the second non- analytic term in the energy expansion of 
the transport scattering time. 

From Eq. (j9j) and Eq. jTO]) it follows 



1 



f 



e F 



n(e) T (e F ) e F 



[B{e F ) + 2C{e F )Q{e-e F )}, (11) 



where B(e F ) = [Bi(e F ) + B2] is the sum of two terms. 
One of the terms in square brackets describes the contri- 
bution of form factors and is equal to 



Eqs. (|16ll7p show us that within the accuracy of our 
calculations the diffusion thermoelectric coefficients in 
the regime of strong electron-electron interaction are sig- 
nificantly lower than the values expected for the case of 
a normal metal with weak interaction. Similar behavior 
is observed in experiment. — Moreover, there is a contri- 
bution which leads to the sign change of the thermoelec- 
tric coefficient L^ 2 at very low electron densities, reach- 
able in remote doped Si/SiGe heterostructures, where the 
Coulomb interaction is strong. 

In our model the thermopower at low temperatures is 
given by the following expression 



Bi{e F ) = £(15-33- 



ei(2k F )' 



(12) 



Sd 



n 2 k 2 T 

■—-^-[l-B(e F )-C(ep)},k B T<^ep. (18) 
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where b = (48irm±e 2 N* /e^) 1 / 3 , TV* = N dep i + ±±n , and 
!f < 1. 

The second term B2 is caused by the local-field cor- 
rection and is independent of electron density and equal 
to 



Bo 



1 



5 3 / 2 g v -5' 



(13) 



For g v = 2 we have B2 — 0.06 and for g v = 1 we have 
B2 — 0.16. The function C(e F ) is expressed by: 



C{e F ) ~ 1 



1 



ei(2k F )' 



(14) 



This C(e F ) term describes the non- analytical behavior 
of the energy expansion of the transport scattering time 
and will be the origin of anomalous temperature depen- 
dencies. The transport scattering time is then given by 



T Q {e F )(l 



e F 



£f 



[B(e F ) + 2C(e F )0{e-e F )}) 

(15) 

After substitution of Eq. (|15l) into Eq. ([2]) we obtain 
the final result for g v — 2 



rl2 



-0.06- 



15 - 33^ k F ireT (e F ) 2 



ei(2k F ) 



k%T. 



It clearly shows the importance of the Coulomb interac- 
tion by the factor C{e F ) containing explicitly and in ana- 
lytical form the interaction potential including exchange 
and correlation. The appearance of this factor C{e F ) in 
the formula for the thermoelectric power of an interacting 
2DEG is the new result of our paper. This factor does 
not appear in three-dimensional systems because there 
the non-analytical term in the scattering time is much 
smaller. But we expect that a similar term should exist 
in the one-dimensional interacting electron gas. 
The parameter B{e F ) also can be written as 



(16) 



Note that the "ideal electron system" contribution and 
the main term in the non-analytical energy dependence 
of the transport scattering time cancel each other. In the 
limit of very small electron density Eq. (1161) reads 
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' m b \ 2 
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0.45 
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ireT (s F ) 



kiT. 



(17) 



Here effective mass m is renormalized by interaction ef- 
fects and m& is the value for the band mass. 



B(e F ) 



e F ,dT (e) 



To(e F ) 



de 



(19) 



We stress that the derivative must be taken at an energy 
from below the Fermi energy e = Ep in order to avoid the 
non-analytic behavior of the transport scattering time. 
This term B{e F ) is in agreement with earlier results on 
the thermopower.— 

For (100)Si MOSFET systems with intermediate elec- 
tron density (3 — 10) * 10 n cm -2 we evaluate the ther- 
mopower as 



S d 



0.47- 



0.06]. (20) 



It is interesting to compare Eqf20] and Eqf4j At low 
density (r s ^> 1) the absolute value of thermopower is 
strongly decreased and the dependence on the electron 
density is expected to be rather weak, because Sd °c r s . 

Until now we have discussed only the low temperature 
behavior of the thermoelectric coefficients. In experi- 
ment the temperature interval for the diffusion regime 
conventionally is defined as the interval in which the 
thermopower is proportional to the temperature. For 
the thermopower in the diffusive regime we get linear 
corrections, which are anomalous due to the anoma- 
lous screening in two dimensions. We conclude that 
Sd/T = ci +C2T. c\ was calculated in the present paper. 

For more realistic calculations and the comparison with 
real experimental results one needs all disorder parame- 
ters such as the density of impurities, the parameters 
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for the interface-roughness scattering and a realistic form 
of confinement potential (depletion density), 24 Moreover, 
for the comparison between theory and experiment it 
would be useful to measure for the same sample as func- 
tion of density and temperature the conductivity and the 
thermopower. 

The main message of our calculation consist in the re- 
sult that the anomalous screening in two dimensions has 
a much larger impact for the Seebeck and the Peltier ef- 
fect than on the conductivity. The effect is such large 
that at low carrier density all calculations, not taking 
into account this contribution, are unable to predict cor- 
rectly even the order of magnitude for the Seebeck and 
the Peltier effect. The corresponding contributions to 



the conductivity [in in Eq.©] and the thermoconduc- 
tivity [L22 in Eq.([2])] are only small corrections, linear in 
temperature , 17 i 18 and disappear at very low temperature. 

Our theory generalizes earlier theoretical results ob- 
tained for the diffusion thermoelectric coefficients in the 
interacting two-dimensional electron gas. The addi- 
tional contribution we obtain depends strongly on the 
strength of the Coulomb interaction including exchange- 
correlation effects via the local-field correction. 
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Shashkin and V. F. Gantmakher. This work was sup- 
ported by the RFBR, the Russian Ministry of Sciences, 
and by the Programmes of RAS. 
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